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Abstract: We present an update on the study of topological structure of QCD. Issues addressed include a 
comparison between the plaquette and the geometric methods of calculating the topological density. We show 
that the improved gauge action based on blocking transformation suppresses the formation of topologically 
charged dislocations with low action. Using a cooling method we identify the instantons' location, estimate their 
size and density, and calculate the renormalization constant Zq for the plaquette method. 



Issues related to the discretization of the con- 
tinuum expression for the topological density, 
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depend on which of several methods, of which 
we consider the plaquette [ |l|, and geometric 
(Liischer's [ ||), is used to calculate q{x), and 
on the lattice action. The geometric method 
while preserving an integer result for the topo- 
logical charge Q = q{x) gives rise to a diver- 
gent topological susceptibility xt in the contin- 
uum limit due to lattice artifacts called disloca- 
tions [ ||]. We show that these dislocations are 
much better controlled by an improved action 
(lA) derived by approximating the renormaliza- 
tion trajectory(RT) of the (a/S) block transfor- 
mation We also show that cooling [|^ is use- 
ful for quantifying instantons' locations and sizes, 
and for estimating the large multiplicative renor- 
malization Zq{P = 6) for the plaquette method[ 

Our calculations are performed on the existing 
ensemble of 34 [Kr" x 40) lattices at /3 = 6.0 using 
the Wilson action [ M]. For the RGT improved 
action, we use 55 (18 x 36) lattices at /3 « 6.0. 

1. Improved Action: Geometric Method 

On the lattice the action of an instanton goes 
to zero as the core size p ^ 0, rather than 
show the expected logarithmic divergence. The 
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charge, measured using the geometric method, 
jumps from to 1 at some Pc where the lattice 
action Smm is still much smaller than that for 
the classical instanton. These boundary configu- 
rations with Q = ±1 are called dislocations, and 
are lattice artifacts because their size p < a for 
all g, i.e. pc does not scale as expected. In the 
case of the Wilson action (WA) the Boltzmann 
suppression e~^^'^^" is overwhelmed by the en- 
tropy factor e^^^'" and dislocations dominate 
the calculation of xt in the continuum limit. To 
raise Smin we investigate an lA consisting of the 
plaquette in the fundamental, 8 and 6 representa- 
tions and the 1x2 planar loop with couplings [ |^] 

^^.0M:^= -0.12; ^ = -0.12. (2) 
Kp Kp Kp 

In our normalization, the bare coupling is given 
by Kp — C(ix 1)6/(7^ where C(ixi) is the weight 
associated with the plaquette. We estimate 5'^„ 
by the following method. We constrain a cen- 
tral plaquette to be —1 in a SU{2) projection, a 
characteristic of most dislocations that have been 
considered, and then vary the links around the 
plaquette in a collective fashion until a minimum 
of the action is found with a topological charge 
of 1. For such configurations S^^^ = 18.9, even 
larger than the value 87r^/6 for a classical instan- 
ton. A comparison of S^^^ with other actions 
studied in [ ^ is made in Table along with 
C(ixi) for each action. We conjecture that dislo- 
cations are suppressed for improved actions with 
large C(ixi) because the central plaquette near 
the cut locus contributes a significant fraction 
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Table 1 



Minimum action of configurations with Q = 1 



Action 


C(lxl) 


SU(3) Smin 


Ref. 


Wilson 


1.0 


4.5 


[ 


3 




Symanzik 


1.67 


6.1 


[ 


3 




Wilson Impr. 


4.376 


10.7 


[ 


3 




VS RGT 


9.09 


18.9 





to the total action of the dislocation and a large 
C(ixi) raises the action of this central plaquette. 

The results for xt by the geometric method at 
Kf = 10.58 and for the WA [ § at /? = 6.0 are 
given in Table ||. A rough estimate of the scale 
for the lA is a^^ ~ 2.3 GeV, i.e. the same as 
for the WA at /3 = 6.0. xt with the WA lies 
significantly above that predicted by the Witten- 
Veneziano formula[ ^, however, due to the un- 
certainty in scale for the lA we cannot yet isolate 
the role of dislocations by comparing Wilson and 
improved action results or give a reliable estimate 
for xi^- However, the two advantages of the lA 
in the study of topology are (i) raising Smin to 
suppress dislocations and (ii) the lA plaquette ex- 
pectation value (□) — 1.871, is much higher than 
the WA value 1.782, consequently the integration 
routines for calculating q{x) converge about 70% 
faster as is expected of smoother fields. 

2. Identifying instantons via Cooling 

In order to correlate hadronic structure with 
topology we need to identify the shapes and sizes 
of all instantons on the uncooled lattices. Ultra- 
violet (UV) noise in q{x) on the uncooled lattices 
has thwarted our efforts to devise a method to do 
this. We therefore adiabatically[ ^ cool a subset 
of lattices (32 WA lattices and 28 lA lattices) for a 
total cooling time r = 10 and measure q{x) using 
the plaquette and Liischer's method for cooling 
times r = 0, 0.3, 1, 3, 6, 10. The results are shown 
m Table |, using Zq{t) as estimated below. Note 

Table 2 



Topological susceptibility by geometric method 



Action 


1/4 

Xt a 


Xt (MeV) 


Ref. 


Wilson 
V3 RGT 


0.114(8) 
0.0975(46) 


262(18) 
224(11) 


(this work) 




10 20 30 



Figure 1. Total Number of charges for lA with 
<lmin{T — 10) — 1 vcrsus the effective radius "r^jf- 

that for r > 1 the two methods give consistent re- 
sults. The complete lack of dependence of xT°"^ 
on T is probably fortuitous. 

On the T = 10 cooled lattices we identify hyper- 
spherical volumes of radius < r^^^. containing 
a total net geometric charge of magnitude greater 
than qmin (we generally set qmin = 1) and con- 
sider these volumes to contain an instanton or 
anti-instanton. Our search routine first isolates 
small instantons and excludes their volumes from 
future search. As a result, if there are overlapping 
instantons, larger instantons {rmax ~ 30) lose 
the spherical shape we normally associate with 
classical instantons, therefore for each instanton 
volume containing v points we define an effective 
radius r^ff = (2v / n'^Y ^ ^ . Figure |l| shows the 
distribution of instanton volumes for lA in terms 
of with r'^^x = 35. Most are in the range 
4a < Teff < 5.5a, or about 0.35 to 0.5 fm. Tak- 
ing the mean size to be re// = 5 and assuming 
a uniform distribution, we estimate the core size 
of classical instantons to be p ^ 0.3 fm. On WA 
lattices we find a total of A^_|_ = 180 instantons 
and 7V_ = 200 anti- instantons, corresponding to 
an instanton density of 1.3fm^^. Similarly, for 
the I A we find N+ = 440, iV_ = 421, which gives 
a factor of two larger density of 2.7 fm~^. 

The average charge {q)y{T) in volumes identi- 
fied with instantons at t = 10 is shown in Fig. ^ 
versus the cooling time. For the lA action we 
find that ((?)^^(t) decreases monotonically with r 
which we interpret as (pair) annihilation of small 
instantons during cooling. The behavior with 
the WA at small r is not clear, but for r > 0.3 
{q)^^{T) shows a similar decrease. 



We use the correlation between the average 
charge inside the instanton volume as measured 
by the plaquette and geometric methods to de- 
fine the renormalization constant to be Zq{t) = 
(9)S''"^-(T)/(g)f°™ (T = 10). For uncooled lat- 
tices we find Z^^(/3 = 6.0) = 0.158(13), and 
Z^'^(/3 R:! 6.0) = 0.230(9). On varying the 
search parameters between 0.8 < qmin < 1.0 and 
24 < r^ax ^ 35, the Zq ranges roughly between 
0.14 and 0.17 (WA), and 0.21 to 0.24 (lA), i.e. it 
has a rather weak dependence on the definition of 
instanton volumes. Including this dependence as 
a systematic uncertainty, we estimate the renor- 
malization constants to be 



Wilson 



Z^^ = 0.16(2), 



Zlf- = 0.22(2). 



(3) 



Our estimate of Z'q^ is consistent with the value 
0.18(2) obtained by Alles et al. [ ^ using a heat- 
ing method. The larger value of Zq^ is consis- 
tent with the expectation that the lA is a better 
approximation of the continuum theory (where 
Zq = 1). We use this estimate of Zq{t) to com- 
pile the data presented in Table ^. 

Our cooling procedure does not guarantee that 
all physical instantons survive. We assume that 
instantons that survive are physical and their 
properties, as a function of cooling time, can 
be studied once their location is known. The 
data show that instantons identified on cooled 
lattices can be tracked back to the uncooled lat- 
tices, where UV noise prevented their identifica- 
tion. Systematic errors due to the uncertainty in 
the scale and the non-uniqueness of the search al- 
gorithm have not been fully included in the rough 
estimates of p, density, and Zq(t) given above. 

We are refining this analysis in several ways. 
We are calculating the string tension for the lA 
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Improved action xj'^ (MeV) (28 lattices) 
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Figure 2. Average charge {q)v{T) contained 
within instanton volume versus r for geometric 
(crosses) and plaquette (diamonds) methods. 

lattices to get a reliable estimate of the scale. We 
are redoing the lA analysis with r = 30 as the 
starting point because the data in Fig. ^ show 
that {q}v{T = 10) with the geometric method 
has not stabilized. We are investigating alter- 
nate schemes for identifying instanton volumes, 
and are generating quark propagators to corre- 
late hadronic properties with instantons for the 
lA configurations. Once the signal in all these 
quantities is under control we plan to do the cal- 
culation at weaker coupling to check scaling. 
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